Periodic orbit theory and spectral rigidity in pseudointegrable systems 
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We calculate numerically the periodic orbits of pseudointegrable systems of low genus numbers g that arise 
from rectangular systems with one or two salient corners. From the periodic orbits, we calculate the spectral 
rigidity As(L) using semiclassical quantum mechanics with L reaching up to quite large values. We find that 
the diagonal approximation is applicable when averaging over a suitable energy interval. Comparing systems of 
various shapes we find that our results agree well with A3 calculated directly from the eigenvalues by spectral 
statistics. Therefore, additional terms as e.g. diffraction terms seem to be small in the case of the systems inves- 
tigated in this work. By reducing the size of the corners, the spectral statistics of our pseudointegrable systems 
approaches the one of an integrable system, whereas very large differences between integrable and pseudointe- 
grable systems occur, when the salient corners are large. Both types of behavior can be well understood by the 
properties of the periodic orbits in the system. 

PACS numbers: 05 .45. -a, 



I. INTRODUCTION 



The motion of a classical particle in a billiard system can 
show regular, chaotic or intermediate behavior, depending on 
the billiard geometry. In a chaotic billiard, the motion is er- 
godically extended over the whole energy surface in phase 
space and two particles whose trajectories are very close at 
the beginning, diverge exponentially from each other. If the 
system is integrable on the other hand, the motion of the bil- 
liard particle is restricted to a two-dimensional torus in phase 
space and neighboring trajectories diverge only linearly from 
each other. Examples for chaotic billiards are e.g. the Sinai or 
the stadium billiard, whereas rectangular or circular billiards 
are integrable. Between these two limiting cases, there are 
several classes of intermediate systems. 

A potential well of the same geometry as the corresponding 
classical billiard - a quantum billiard - reflects this regular, 
chaotic or intermediate behavior in the statistics and the dy- 
namics of its eigenvalues. The statistics of the eigenvalues 
investigates the static correlations of the eigenvalues and the 
distribution of the distances between consecutive values [1]. It 
can be determined directly by calculating first the eigenvalues 
and then their distribution and the correlations between them. 
It can also be determined at least approximately by calculating 
the periodic orbits of the system and applying semiclassical 
quantum mechanics. 

In this paper, we want to use the periodic orbit theory. 
We focus on pseudointegrable systems [2-4], which are an 
interesting example of an intermediate class. Like in inte- 
grable systems, the motion of a classical particle in a pseu- 
dointegrable system is restricted to a two-dimensional sur- 
face in phase space. However, these surfaces do not have 
the shapes of tori but are more complicated objects with more 
than one hole ("multihandeled spheres"). Examples for pseu- 
dointegrable systems are polygons with only rational angles 
nin/nii, with m,mi G N and at least one m > 1. They are 



classified by their genus number 
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which is equal to the number of holes in the multihandeled 
sphere in phase space. Here, J is the number of angles and 
M is the least common multiple of the m,. The reason why 
those systems are not completely integrable is their property 
of beam splitting. At some points in their geometry, neighbor- 
ing trajectories of particles can be split into two opposite di- 
rections. Fig. 1 shows examples of pseudointegrable billiards 
with genus numbers g = 2 and 3. The beam splitting property 
is demonstrated in Fig. 1(a) at one of the salient corners. 

The first possibility to calculate the spectral statistics starts 
with the eigenvalues. Here, one first looks at the distribu- 
tion P(s) of the normalized distances s, : = (£j+i — £$)/ (s), 
(s) = 1, between two consecutive energy levels e i+1 and 
with the mean distance (s), which has two limiting cases. 
When the systems are integrable, P(s) follows the Poisson 
distribution, whereas the Sj of chaotic systems are Wigner- 
distributed. The distribution P(s) of pseudointegrable sys- 
tems has been found to be intermediate between both [5-10] 
and it is assumed that with increasing g they come closer to 
the behavior of chaotic systems. However, it is not yet clear, 
in which way other system details interfere. 

As a measure for the correlations between the eigenvalues, 
we consider the spectral rigidity As(L) [11], which describes 
the correlations in a normalized energy interval of length L. 
L gives the approximate number of energy levels in the con- 
sidered interval and will in the following be called the argu- 
ment of A3. We start from the integrated density of states 
N(e) = ^2n =1 B(e — e„) of the normalized ("unfolded") en- 
ergy levels, which is a staircase and can be approximated by a 
straight line. A 3 (L) is defined as the least square deviation, 
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where Min ri r2 means that the parameters n and r 2 are cho- 
sen such that the line n + r 2 e is the best fit of N(e). The 
average (. . .) is an energy average, carried out over many 
different values of Eq in an interval AE. AE should not 
be confused with the argument L, which gives the length of 
the considered energy interval. The limiting curves for not 
too large L are A 3 (L) = L/15 for integrable systems and 
A 3 (L) = ln(£)/7r 2 - 0.07/tt 2 + Oil- 1 ) for the ensemble 
of Gaussian orthogonal matrices (GOE) [1,12], which serves 
as a generally accepted good limit for chaotic systems. For 
pseudointegrable billiards, A 3 (L) has been found intermedi- 
ate between both (see above). 
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FIG. 1.: Shapes of the pseudointegrable systems considered in this 
paper: (a) the one-step system (Xi, Yi; X2, Y2) with genus number 
g — 2 and (b) the two-step system (Xl, Yi; X2, Y2; X3, Y3) with 
g = 3. In (a), the beam splitting of two initially neighboring tra- 
jectories at a salient corner is demonstrated and in (b) trajectories of 
three "neutral" orbits (between two parallel walls) are indicated by 
arrows. For the lengths Xi, Yi considered in this paper see Tab. I. 

A second possibility to calculate A 3 (L) is given by the pe- 
riodic orbit theory. In pseudointegrable systems, all periodic 
orbits form families of equal lengths and the starting point of 
an orbit can always be shifted to at least one direction along 
the boundary without leaving the family. The simplest and 
shortest orbit families are the "neutral orbits" [13] that bounce 
between two parallel walls. Using semiclassical quantum me- 
chanics, setting h = 2m = 1 and neglecting additional con- 
tributions coming e.g. from diffractive orbits, A 3 (L) under 
Neumann boundary conditions is given by [12,14-16], 



occupied by the respecting orbit families and the function 
Hij = F{ Vl - Vj ) - F{ yi )F{y ) - 3F'( yi )F'{ yj ), where 
F(y) = (smy)/ y and primes denote differentiation. The ar- 
gument L enters via yi = (Lti) / [A^/Eq {d)) with the mean 
level spacing 1/ (d) = At:/ A and the system area A. As in 
Eq. (2), the average (. . .) in Eq. (3) is carried out over differ- 
ent energies Eq in an energy interval of width AE. 

In this paper, we basically want to use Eq. (3) to calcu- 
late A 3 (L) and see in which way the different system details 
apart from the genus number g influence its behavior. It can 
be seen by Eq. (3) that the number of orbits in the different 
length intervals [£, £ + A£] as well as the corresponding areas 
a(£ ) are the important quantities to investigate this question. 
We therefore carefully calculate these quantities for systems 
with different lengths and widths (Xi, Yi) of the different seg- 
ments (cf. Fig. 1) and see how the behavior of A 3 changes by 
varying the systems. 

The behavior for large L is determined by the short orbits 
and we will sometimes use the neutral orbits for crude approx- 
imations. By choosing the segments X2 — > X\, Y 2 — * Y 2 , 
the orbits approach the ones of a rectangular system, which 
is integrable. For large differences between X\ and X2 and 
between Y\ and Y 2 , on the other hand, the orbit families can 
become very different from those of the rectangle. 

The paper is organized as follows: In section II, we ex- 
plain how the lengths and areas of the periodic orbits are cal- 
culated and show the results. In section III, we discuss the 
applicability of the diagonal approximation, where only terms 
of li — £j are considered in Eq. (3). In section IV we finally 
show the periodic orbits results for A 3 (L) for many different 
systems and compare them to the eigenvalue statistics. Some 
of the considered systems are very close to integrability, while 
other systems possess very pronounced salient corners and we 
show how A 3 (L) is influenced by these system details. Fi- 
nally, in the conclusion section V, the influence of neglected 
terms as e.g. the diffraction terms and possible deviations be- 
tween the eigenvalue statistics and the periodic orbits results 
are discussed. 



II. CALCULATION OF THE PERIODIC ORBITS 

It has been shown in [14,15] that for large £, the prolifera- 
tion rate N(£), i.e. the number of orbits with lengths smaller 
than £ grows quadratically with £, 



N{£) = irb £ 2 /(a(£)}, 



(4) 
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The double sum is carried out over all orbit families (in- 
cluding repetitions, but only in forward direction) of lengths 
£i and £j, a, and aj are the areas in phase space that are 



where (a(£)) = J2n <e a i/J2n <£ 1 i s tri e average area in 
phase space, occupied by the orbits with lengths smaller than 
£ and b is a constant, depending slightly on the details of the 
system, 



b S3 (l/2n)J28(£- £i) ai /£i 
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(the sum going over all periodic orbits including repetitions). 
It has been found that b = 1/4 for integrable systems 
and slightly larger for pseudointegrable systems, whereas 
(a(i)} ~ 4:A in integrable systems and considerably smaller 
for pseudointegrable systems [15]. 
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FIG. 2.: Test of the sum rule (see Eq. 6) for the geometries of Fig. 1 
with different values of Xi,Yj. S(£) is plotted versus I (giving 
straight lines) and bo is calculated from the slopes. The different 
symbols refer to the rectangular system (triangular symbols) and to 
pseudointegrable systems of g = 2 (circles) and g = 3 (squares). 
Open and closed symbols refer to different step sizes (for details of 
the systems see Tab. I). For the rectangle, we find bo = 1/4 and 
for the pseudointegrable systems, bo is slightly increasing with g and 
with the step sizes. 

We have calculated the periodic orbits of our pseudointe- 
grable systems by two different methods. The first method 
uses the fact that one member of each orbit family must ei- 
ther start at a salient corner or pass close to it. Therefore, 
we start at the salient corners and vary the reflection an- 
gle ip between the trajectory and the billiard wall in small 
steps Aip until the orbit nearly closes. Finally we use an 
iteration method [17] to find the exact reflection angle ip. 
The second method calculates ip and the lengths £ v of all 



hypothetical orbits by [18] tan</? 
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U = 2[(E ! ^) 2 + (E j m 3 X J ) 2 
are positive integers and Xi, Yj are the segment lengths as 
shown in Fig. 1. In pseudointegrable systems, due to the 
shielding role of the corners, not all hypothetical orbits re- 
ally occur in each system. Therefore, we have to check which 
trajectories actually return to their starting point within the 
correct length £ v of the trajectory. Both methods basically 
lead to the same results. The agreement, e.g. for the one-step 
system is larger then 97% for the orbits, being reflected at the 
boundaries up to 50 times. Going to more reflections or to sys- 
tems with more steps, the iteration method misses more orbits, 
depending on the value of Aip and the numerical tolerances, 
which makes the second method more reliable and therefore 
better suited for this work. Moreover, the second method is 
faster for our systems of small genus numbers, but we assume 



that the iteration procedure could be favorable in more com- 
plicated systems, where the number of hypothetical orbits can 
become very large. 

Naturally, these numerical procedures are restricted to a 
maximum orbit length, in our case to the first 3000 orbits (for 
tests 12000), including repetitions. This is legitimate as the 
contributions of the diagonal elements in Eq. (3) decay with £ 3 
(the contributions of non-diagonal elements can be neglected 
as will be shown in the next section) whereas the number of 
orbits per length interval A£ increases only with 2£A£ (see 
Eq. (4)). With Hij staying finite, the alternating sum of Eq. (3) 
is convergent. 
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FIG. 3.: Normalized average area (a(£)) in phase space for the pe- 
riodic orbits with lengths smaller than I for the same systems (and 
same symbols) as in Fig. 2. It can be seen that systems of the same 
genus number g form groups of similar (a(£)). For large values of £, 
the areas are saturating. 

First, we verify that our periodic orbits fulfill Eq. (4) and 
investigate how bo depends on the system details. bo can ap- 
proximately be calculated from the sum rule [15] 



S(£) = S(£ ini J = £ ai/ii « 2nb £, 



(6) 



which can be easily verified by replacing the sum by an inte- 
gral over d£, inserting the density dN(£)/d£ with N(£) from 
Eq. (4) and replacing the areas aj by their mean value (a). In 
Fig. 2, we plot S(£) versus I for the pseudointegrable systems 
of Fig. 1 with different step sizes as well as for the rectangular 
system. In all cases, S{£) increases linearly with the upper 
orbit length £ and bo is determined from the slopes by a least 
square fit. We can see that 6 increases slightly with the genus 
number g and the step sizes. 

Next, we calculate the average area (a(£)) in phase space 
for the periodic orbits with lengths smaller than £. In Fig. 3, 
we plot the normalized average area (a(£)} / {A A) versus £. In 
agreement with previous calculations [15,18], the values of 
(a(£)} are saturating for large values of £ and are considerably 
smaller than the value of (a(£)) w AA of integrable systems. 
Moreover, it can be seen that systems of the same genus num- 
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ber g form groups of similar (a(£)}. This can be most easily 
understood by comparing the neutral orbits of Fig. 1(b): If 
we start with a rectangle and disturb it by one salient corner, 
each neutral orbit will split into two different families, each 
of them covering a smaller area than before. The same effect 
also occurs for more complex orbits and is repeated with each 
additional salient corner. So, the average area (a(£)} will be 
reduced with g, whereas the number of different orbit families 
increases. 

Now, we can calculate the proliferation rate N(£) for our 
systems according to Eq. (4) by using the values of 6 an d 
(a(£)) from Figs. 2 and 3, respectively. In Fig. 4 we com- 
pare the values from Eq. (4) (straight lines) to the numbers 
of N(£), obtained by counting the different orbits (symbols). 
The agreement between both curves is excellent. Clearly, 
N(£) increases much faster for higher values of g than for 
lower values (see above). Therefore, also N(£) forms groups 
of very close-lying curves that correspond to systems with the 
same g. Like (a(£)) also N(£) is basically determined by g 
and changes only very slightly by further details of the con- 
sidered systems. 
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FIG. 4.: The proliferation rate N(£) is plotted versus I for five dif- 
ferent systems and two types of calculations. The lines are calculated 
according to Eq. (4), using the values of 60 and (a(tj) from Figs. 2 
and 3, respectively. The symbols (same symbols as in Fig. 2) show 
the numerical data, gained by a simple summation of the different 
periodic orbits (including repetitions). Both curves agree very well. 



III. THE DIAGONAL APPROXIMATION 



As Eq. (3) involves a double sum over £i and £j, it is very 
convenient to use the diagonal approximation, where only 
terms with L = are taken into account. It has been shown 



in [12] that this is justified for integrable systems. The reason 
is that only orbit pairs where i t — £j ~ l/\/AE survive the 
average over an energy interval AE, while for larger phase 
differences, the different terms cancel by the oscillations of 
the cosine function. With N{£) <~ £ 2 and for not too small 
AE, the number of these pairs grows more slowly than their 



contributions decrease and non-diagonal terms can therefore 
be neglected [19]. 

As pseudointegrable systems obey the same kind of 
quadratic proliferation rule as integrable systems (cf. Eq. (4)), 
one can assume that the diagonal approximation should ap- 
ply by the same reasons. However, by calculating the related 
quantity 
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doubts on the validity of the diagonal approximation have 
been expressed [15], because due to the different proportion- 
ality factors, N(£) nevertheless increases faster than in inte- 
grable systems. Deviations between the diagonal approxima- 
tion and the full summation have been found in [15] by inves- 
tigating numerically the integral 
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where £ max = £/ (A-K\fE~o(d)) and the energy average runs 
again over several values of E . 
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FIG. 5.: 1(1) for a rectangle (upper three curves) and a two-step 
system (lower three curves) (for details of the systems see Tab. I). 
In both cases, the diagonal approximation 1(1) is close to a straight 
line (solid line), whereas the full sum oscillates around it (dotted and 
dashed lines). The oscillations are significantly reduced by averag- 
ing over more energy values. The sums are carried out over the first 
3000 orbits, including repetitions. 

Here, we want to take a closer look to these deviations and 
therefore investigated I(£) for our systems by calculating ex- 
pression (8) first in diagonal approximation and second by 
carrying out the full sum. When some orbit lenghts are de- 
generated, we include all pairs of orbits with £i = £j into the 
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diagonal approximation. This is equivalent to combining or- 
bits with ii — lj to one single orbit family with area a, + aj . 

The results are shown in Fig. 5. The upper three curves 
show the rectangular system and the lower ones the two-step 
system. Both systems are evaluated (i) in diagonal approxi- 
mation (solid lines) and with the full sum and an average over 
(ii) 160 (dashed lines) and (iii) 20 values of Eq (dotted lines), 
showing small and large fluctuations around the solid line, re- 
spectively. The results for the one-step system, that are qual- 
itatively the same, are omitted for a better overview of the 
curves. 

Contrary to [15], we do not find a crossover value for t in 
pseudointegrable systems above which the diagonal approx- 
imation breaks down. Instead we find fluctuations for I(£) 
around the diagonal approximation In (I) that depend on the 
energy average and increase with £. They can become quite 
large and seemingly distant from Id{£)- However, these fluc- 
tuations occur in both, the integrable and the pseudointegrable 
systems and are significantly reduced by the average proce- 
dure. Therefore, the diagonal approximation seems to be valid 
for both kinds of systems, as long as the number of different 
E values as well as their interval AE is not too small. 

IV. SPECTRAL RIGIDITY 

We finally investigate the spectral rigidity A 3 (L). For in- 
tegrable systems, it is known that A^(L) shows linear behav- 
ior for small L (see section I) and reaches a plateau for large 
L, the height of the plateau being determined by the smallest 
periodic orbit of the system. We now use the periodic orbit 
theory to find A^(L) for our pseudointegrable systems. 

One may ask, if the grouping of (a) and N(£), observed in 
section II also leads to an arrangement of the A3 -curves ac- 
cording to their genus numbers. This can be roughly estimated 
from Eq. (3) by again investigating the neutral orbits. By go- 
ing from the rectangle to the one-step system, the neutral orbit 
families splits into two new ones. With a characteristic step 
size of half the system size, one of the two new orbits occu- 
pies half the area ai than before while the area occupied 
by the other one is even smaller. Considering the length of the 
new orbits as roughly constant (which is of course only true 
for the first one), their common diagonal contribution to A3 
is oc (a\ + a\), roughly half the contribution oc (2ai) 2 of the 
respective rectangular orbit. A similar discussion applies also 
for the more complicated orbits and for higher genus numbers. 
Therefore, A 3 of a rectangle should be roughly diminished by 
a factor of two, when introducing a step of half the system size 
(g = 2), and by a factor of three, when introducing two steps 
of equal size (g — 3). 

However, the situation is very different when the step size 
is small compared to the system size. In this case, the contri- 
bution ai of the first new orbit is very close to the contribution 
of the rectangular orbit, whereas the contribution becomes 
negligible. Therefore, for small step sizes, the splitting of one 
orbit family into two new ones can be neglected and we expect 



A3 curves very close to the one of the rectangular system. In 
this case, the A 3 curves should not be ordered according to 
their genus numbers. 

We now want to investigate these assumptions by our cal- 
culations. First, A 3 (L) is calculated by the diagonal approx- 
imation (open symbols) and by the full sum (solid symbols) 
and is plotted in Fig. 6 versus L for the rectangular system 
(squares) and for our systems with g = 2 (circles) and .9 = 3 
(triangles). The energy average was carried out over 160 val- 
ues of Eq and the number of orbits, taken into account, is fixed 
to 3000 including repetitions. We carefully checked that the 
inclusion of even more orbits does not change the results of 
any of the considered systems significantly. The results can 
be summarized as follows: 

The calculated A 3 (i) values for the rectangle are in quite 
good agreement with the theoretical slope of L/15 for small 
L. Also in the pseudointegrable cases, A 3 (i) seems to in- 
crease linearly for small L and saturates for larger L, but the 
slopes are smaller and decrease with increasing genus num- 
bers, thereby coming closer to the A 3 curve of the chaotic 
systems (lowest curve). The agreement between the diagonal 
approximation and the full sum is excellent. 

12 i , , , , , , , , 1 
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FIG. 6.: As(L) is calculated by the diagonal approximation (open 
symbols) and by the full sum (solid symbols) and plotted versus L 
for different systems of genus number g — 1 (squares), g — 2 (cir- 
cles) and g — 3 (triangles) (for details see Tab. I). The theoretical 
behavior for integrable and chaotic systems is indicated by the upper 
and the lower line, respectively (dotted lines without symbols). 

Next, we investigate the dependence on g. When the step 
sizes are roughly equal in size, the curves of Fig. 6 show a sys- 
tematic behavior on g, with A 3 (i) decaying towards smaller 
values with increasing g. The curves for the rectangle are 
highest and most curves for the g = 2 systems are above those 
with g = 3. As expected from the discussion above, a system- 
atic dependance on g occurs also for the lowest curves of each 
genus number (refering to steps that split the system into parts 
of roughly equal size). These curves lie roughly by a factor 
of 2 (for g = 2) and of 3 (for g — 3) below the curves of 
the rectangle, which could be surprisingly well estimated by 
the behavior of the neutral orbits (see above). However, for 
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very small step sizes, this behavior changes and only slight 
deviations of A 3 from the curve of the rectangle occur. This 
can be seen on the upper three curves, which refer to g = 1, 2 
and 3, but to very small step sizes and means that for given 
g all curves between the Lj 15 line and a minimum curve can 
be found by changing the step sizes appropriately (the slope 
of the minimum curve depending on g). Therefore, beside 
the genus number g, the details of the specific system play an 
important role as well. 

It is however interesting to note that with fixed g and in- 
creasing step sizes, the curves decay quite rapidly towards 
their minimum curve and that for small steps sizes, the 
changes in size have the strongest effect. For the top four 
curves of the one-step systems (circles), for example, the step 
size was varied over 10% of the system size. The decrease in 
the slope is over 30%, and more then half of the maximum 
decrease. The very large change in step size from the 4th to 
the 5th curve, on the other hand, has a comparatively small ef- 
fect. This is the reason that most of the curves shown in Fig. 6 
show indeed a systematic dependance on g. 

Another interesting quantity is the plateau value: For large 
L, the values of A3 saturate to a plateau, whose height is de- 
termined by the small orbits. In order to estimate this height, 
one normally uses the approximation [12] 

r°° r ]f Tf 

=2w] o ^ mG{ rnn ] (9) 

with $ from Eq. (7), t = £/(2E Q ) and the so-called orbit se- 
lection function G(x). G(x) can roughly be approximated by 
a step function which is equal to 1 in the plateau regime. For 
integrable systems, (f>(£) = (d)/(2n) for I > £ mln (where 
4nin is the length of the shortest orbit of the system) and 
(f)(1) = for t < £ m i n , which makes the calculation of A 3 
relatively easy [12]. In the plateau regime, the integral (9) de- 
pends only on the lower integration limit that is determined by 
the length of the shortest orbit of the system. 

For our pseudointegrable systems, things become more 
complicated. We can find <p{£) from the slope of I(£) in Fig. 5 
in diagonal approximation, but the values differ for the dif- 
ferent systems. This means that the behavior of the plateau 
depends on the length of the shortest orbit as well as on the 
system dependent value of &(£). As both quantities can be 
varied independentely from each other, there is no easy ex- 
pression for the plateau value as in the case of the integrable 
systems. 

Finally, in Fig. 7, we compare Ag°(L) from the peri- 
odic orbit results (open symbols, diagonal approximation) to 
A3 L (L) from the eigenvalue statistics (filled symbols). To 
this end, we have calculated the eigenvalues under Neumann 
boundary conditions by the Lanczos algorithm (which in- 
volves a discretization of the lattice) and obtained A^' h (L) 
by using the technique as derived in [20]. Our results are in 
line with previous numerical calculations on the two-step sys- 
tem, where the eigenvalues have been obtained for small val- 
ues of L by numerical diagonalization algorithms ( [5], [15] 



and [18]) and by the boundary element method ( [6] and [7]). 
For the rectangle, also the exact eigenvalues for a continu- 
ous system are used for comparision. For the rectangle, the 
A3 curve as obtained from the exact eigenvalues (solid dia- 
monds) agrees very well with Ag° from the periodic orbit 
theory. This had to be expected for a system without diffrac- 
tive corners. The A3 L curve, obtained from the numerical 
eigenvalues of a discretized systems lies somewhat higher and 
the deviation gives us the error bar arising from the discretiza- 
tion. 

The results of several pseudointegrable systems of g = 2 
and g = 3 systems are shown in Fig. 7 as well. Also in these 
cases, the agreement between A3 L and A3 is quite good. 
Again, Ag L lies slightly above A3 , which can in princi- 
ple be due either to the discretization or to the neglection of 
higher-order terms, as e.g. the diffraction terms. The devia- 
tion seems even smaller than for the rectangle. However, by 
calculating rectangles of different sizes, we found that the de- 
viations between A3 L and A3 are in most cases smaller than 
the one shown here, which we therefore consider as the upper 
limit of the error due to discretization. 
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FIG. 7.: A^°(L) from the periodic orbit theory in diagonal ap- 
proximation (open symbols) is compared to Aj' L calculated from 
the eigenvalues under Neumann boundary conditions obtained by 
the Lanczos algorithm (solid symbols) for different systems. For the 
rectangle, A| as calculated from the exact eigenvalues for a conti- 
nous rectangluar system is shown as well (solid diamonts). Different 
symbols represent the genus numbers g = 1 (squares), g — 2 (cir- 
cles) and g = 3 (triangles). The theoretical behavior for integrable 
and chaotic systems is indicated by the upper and the lower line, re- 
spectively (dotted lines without symbols). 

On the basis of the present data, it seems that the contri- 
bution of the higher-order terms is quite small and at most in 
the same order of magnitude as the errors due to discretiza- 
tion. It is also interesting to note that the deviations do not 
increase, when going from g = 2 to g = 3, which doubles the 
number of diffractive corners. This may also be a hint that the 
diffraction terms are not very important in our systems, where 
the diffraction terms are arising from the corner scattering. 
The situation seems to be different in systems with point-like 
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scatterers, where the spectral statistic is known to be changed 
drastically [21]. It will be interesting to investigate systems 
with higher g to see, if this assumption holds. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have calculated the spectral rigidity A 3 (L) 
for various pseudointegrable systems by applying periodic or- 
bit theory. By averaging over enough energy values taken 
from a not too small interval, we found that the diagonal ap- 
proximation is applicable. Most A 3 (L) curves decrease with 
increasing genus number g, but other details of the geometry 
play an important role as well. In particular, when the salient 
corners are very small, all curves approach the one of the rect- 
angle, independentely of their genus numbers. 

The behavior of the different A 3 (L) curves can be under- 
stood in terms of the periodic orbit families. Estimating the 
behavior of the orbits from the neutral orbits that bounce be- 
tween two parallel walls gives already a crude approximation 
of the lowest curves of a given g that works surprisingly well. 
This estimation also shows that the lowest curves occur for 
systems with salient corners of roughly equal size. If the cor- 
ners of the systems are very small, on the other hand, the op- 
posite behavior occurs and the A 3 (L) curves lie very close to 
the one for a rectangle. Also this can be quite well understood 
by looking at the periodic orbit families. The ones that pos- 
sess the same lengths as the rectangle occupy in this case the 
dominant areas in phase space, while the orbit families that 
differ from the ones of the rectangle become negligible. 

Finally, we discuss which effects have been neglected and 
could lead to deviations from our results. First, there are var- 
ious attempts in literature to go beyond periodic orbit theory, 
especially the attempt to take the influence of diffractive or- 
bits into account. These orbits start at the singularities, in 
our case at the salient corners, and their contributions to spec- 
tral statistics have been treated in [22]. They are in general 
non-negligible (except for L — > 0), but smaller than the peri- 
odic orbit contributions. Second, there might be an effect of 
the larger orbits beyond the first 3000 ones. However, as we 
discussed in section III, their influence should be negligible. 
Indeed, the difference between A 3 calculated by taking 3000 
or 12000 orbits (including repetitions) was hardly visible. 

We can get hints about the importance of the diffractive or- 
bits by comparing the A 3 curves obtained from the eigenval- 
ues and the ones obtained from the periodic orbit calculations. 
As shown in the last section, the agreement between both is 
quite good. Moreover, the difference between the periodic or- 
bit results and the eigenvalue results is roughly the same in 
all considered systems, independentely of number and size of 
the salient corners. This indicates that the (small) difference 
between both methods is rather due to the discretization of the 
eigenvalue calculations. Therefore, in agreement with [15] we 
think that the influence of diffractive orbits to spectral statis- 
tics in our systems is only small. Of course, it is possible that 



it increases when still more corners and/or additional diffrac- 
tion terms, arising e.g. from flux lines or point like scatterers, 
are present. 

Moreover, it will also be very interesting to investigate the 
connection between the A 3 curves and the periodic orbits for 
more complicated systems with higher genus numbers. Intro- 
ducing more and more corners leads to more and more differ- 
ent orbit families with quite similar lengths, but rather small 
areas in phase space. When many - even small - corners are 
present, it is quite obvious that the original orbits of the rect- 
angle are considerably disturbed, i.e. that their area in phase 
space becomes much smaller. Therefore, contrary to the rela- 
tively simple systems of this work, one should expect a con- 
siderable change of A 3 . Many numerical works have shown 
that the curves are shifted towards the one of chaotic systems 
and it will be very interesting to rely this shifting also to the 
periodic orbit families of the respective systems. 
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TABLE I.: Table of the systems used in this work. X, and Y t are the lengths of the different segments of the systems (also lengths of the 
"neutral orbits") as shown in Fig. 1. The last columns indicate in which figures the different systems occur and the numbers given in these 
columns enumerate the curves from top to bottom. 



8 



